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Abstract. We construct semi-infinite billiard domains which reverse 
the direction of most incoming particles. We prove that almost all parti- 
cles will leave the open billiard domain after a finite number of reflections. 
Moreover, with high probability the exit velocity is exactly opposite to 
the entrance velocity, and the particle's exit point is arbitrarily close to 
its initial position. The method is based on asymptotic analysis of sta- 
tistics of entrance times to a small interval for irrational circle rotations. 
The rcscalcd entrance times have a limiting distribution in a limit when 
the number of iterates tends to infinity and the length of the interval 
vanishes. The proof of the main results follows from the study of related 
limiting distributions and their regularity properties. 



1. Introduction 

The present paper is motivated by the problem of constructing open bil- 
liard domains with exact velocity reversal (EVR), which means that the 
velocity of every incoming particle is reversed when the particle eventu- 
ally leaves the domain. This problem arises in the construction of perfect 
retroreflectors — optical devices that exactly reverse the direction of an inci- 
dent beam of light and preserve the original image. A well-known example 
of a perfect retroreflector is the Eaton lens [3], [12] which is a spherically 
symmetric lens that, unlike our model, also reverses the original image. A 
second application lies in the search for domains that maximize the pressure 
of a flow of particles |10| : for a particle of mass m > 0, which moves towards 
a wall with velocity v, the impulse transmitted to the wall at the moment 
of reflection is equal to 2m|v n |, where v n is the normal component of v. It is 
maximized when v = v n , i.e. when the direction of the particle is reversed. 

We construct a family of domains D e , for which EVR holds up to a set of 
initial condition whose measure tends to zero in the limit e — > 0. 
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The domain D e is the semi-infinite tube [0, oo) x [0, 1] with vertical barriers 
of height e/2 at the points (n, 0) and (n, 1), n G N as illustrated in Fig. 1. 
Inside the domain the particle moves with the constant speed and elastic 
reflections from the boundary. Since the kinetic energy of the particle is 
preserved, we can assume that the speed of the particle is equal to one. 

The motion of the particle is determined by the point y in G [0, 1], where 
it enters the tube and the initial velocity v- m = (cos(inp), sm(inp)) at this 
point. The measure P on the initial conditions {y- m ,tp) considered below 
is a Borel probability measure absolutely continuous with respect to the 
Lebesgue measure on Q = [0, 1] x [—1/2, 1/2]. 



Theorem 1. For every e G (0, 1) there exists a set 12(e) C Q of full Lebesgue 
measure, such that for every (y in ,<p) G 12(e), the particle eventually leaves 
the tube. 

The position and the velocity with which it leaves the tube are denoted 
by (y ut,v out )- By Theorem [TJ for every e G (0,1) the functions y ont = 
y ut(yin,Vin) and v out = v out (y in ,v in ) are defined P-almost everywhere. 

Theorem 2. For any 5 > 0, 

(1.1) F{(y in ,v in ) : v out = -Vm, \y ut - Vm\ < 8} -> 1 as e -> 

Theorem [2] follows from the results on the existence of certain limiting 
distributions for the exit statistics of the billiard particle as e —>■ 0. Below 
we formulate these results as Theorem [3] and Theorem 01 In the last section 
of the paper we show how they imply Theorem [2j 

Let Q e = Q e (yin,Vm) be the number of reflections from the vertical walls 
before the particle leaves the tube. Let T £ = T £ (y in , t> in ) be the time that 
particle spends inside the tube. By Theorem [TJ both Q e and T £ are finite 
P-a.e. 

Consider also a bi-infinite tubular domain similar to the one described 
above. It consists of two horizontal lines at the unit distance from each 
other and a one-periodic configuration of vertical walls of height e/2. 




Figure 1. The Model 
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Let x be the horizontal coordinate and assume that the particle starts at 
x = 0. Let £° = and ^ £ Z be x-coordinate of the particle at the moment 
of fc'th reflection from a vertical wall. Since the tube is now bi-infinite, 
is a discrete time process on Z, defined for any k E N. We also define a 
continuous version of this process: {£ e (t)} is the projection of the trajectory 
of a billiard particle in the bi-infinite tube to the x-axis normalized to have 
constant speed \je. 

Theorem 3. (1) The process {e£g} converges in distribution (w.r.t. F) 
to a stochastic process as e — > 0. 
(2) There exists a limiting probability distribution function G : N — > [0, 1] 
sitc/i i/iai /or ever?/ fc£N, P{Q e (?/m, ^in) = k} —> G{k) as e —* 

The second part of Theorem [3] says that for the limiting stochastic process 
with probability one there exists k G N, such that £ fe < 0. Similar 
results are true for the continuous process {£ e (s)} as well: 

Theorem 4. (1) The process {e£ e (s)} converges in distribution w.r.t. P 
to a stochastic process £(s) as e — > 0. 
(2) There exists a limiting probability distribution function H : IR>o — > 
[0,1], such that for every t > 0, P{eT e (?/i n , t> in ) < £} — >• as 
£ -> 0. 

2. Reduction to Circle Rotations and Point- Wise Exits 

We first reformulate the problem in terms of circle rotations. 
Let us identify [0, 1) with S 1 = R/Z. For a£l, let R a : S 1 -> S 1 be the 
circle rotation by angle a : 

= £ + a mod 1. 

Always assume that a £l\Q. 
Let I £ = [-e/2,e/2] C S l . 

We define several sequences measuring the return times to the interval I ei 
which will be used throughout the proofs. The hitting times = m^(x, a), 
k = 0, 1, 2, . . . are defined for x G S 1 by: 

m° £ = 0, m k £ (x) = min{Z > m k ~ x : R l a x G J £ } 

The sequence n k e = n k (x,a), k = 1,2,... of relative return times to the 
interval I e is defined for x G S 1 by: 

n\ = m k e {x) - m k ^ 1 (x) 

We shall also use the sequence defined in the introduction as the 
sequence of the horizontal coordinates of points of the reflection from the 
vertical walls. 
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Note that if x = y- m , and a = tan(7r<£>), then n\ (x) is the distance between 
horizontal coordinate of the place of the (i — l)'st and the i'th reflections of 
the particle from vertical walls. Therefore, 

g = nl-n 2 £ + ... + (-l) k+1 nt 

Q e = Q £ ( x , a) = min{j G N : n\(x) - n 2 £ (x) + . . . + (-l) j+1 n E (x) < 0} - 1 
Let n k = (n\, n k e ) T ', fh k = (m\, . . . , m k ) T and £ k = , . . . , t; k ) T , then 

(2.2) l k = An k , and m k = Bn k , 

where A and B are two k x k matrices with 



_ J 0, if i < j, _ j 0, if i < j, 

a m - \ (-iy+\ if i > j ' ana ~ \ i, if i > j 

The probability measure P on the initial conditions (y- m , ip in ) G [0, 1] x 
[—1/2,1/2] for the billiard particle induces a probability measure on the 
initial conditions (x, a) G [0, 1) x [0, 1) ~ T 2 for the circle rotation R a , 
which is absolutely continuous w.r.t. to the Lebesgue measure on T 2 and 
which will be also denoted by P. 

We now prove Theorem [l] 

Let T a ^ e : I E — > I £ be the map induced on I £ by the circle rotation R a : 

f a;£ ( X ) = Rf*(x) 

Proposition 5. For every e G (0, 1) there exists a set of full Lebesgue 
measure A(e) C S 1 , such that for every a G A, the map T a ^ e is weakly 
mixing. 

Proof. The proof of Proposition will follow from a combination of results of 
[l\ and [2]. 

For every e > there exists a full Lebesgue measure set A'(e) C S 1 , such 
that for every a G A'(e), the corresponding map T a>£ is an interval exchange 
transformation of three intervals of combinatorial type (3, 2, 1). 

Recall Property P introduced by Boshernitzan in [I]: 

Definition 1. A set A C N is called essential, if for any integer I > 2 there 
exists A > 1, such that the system 

n i+ i > 2ui, for 1 < i < I — 1, 
ni < An 1; 

Hi G A, for 1 < i < I 
has an infinite number of solutions (jii, ri2, ■ ■ ■ n{). 

Let m n (T Qj£ ) be the length of the smallest interval of continuity of T™ £ . 
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Definition 2. An interval exchange map T Qj£ has property P, if for some 
5 > the set 

A(a,e,5) = {neN\m n (f a , £ ) > -} 

is essential. 

Proposition 6. fll], Theorem 9.4 (a)) For every e G (0, 1) there exist a full 
Lebesgue measure set A(e) C S 1 , such that for every a G A, the map T afi 
has property P. 

By Theorem 5.3 of [2J, property P implies weak-mixing for an interval 
exchange of three intervals with combinatorics (3, 2, 1) (and more generally, 
for any combinatorics of a so-called W— type, see [2]). 

This implies Proposition |5j □ 

The next two statements are well-known. We include their proofs to keep 
the exposition self-contained. 

Lemma 7. For every e G (0, 1) there exists a set of full Lebesgue measure 
A(e) C S l , such that for every a G A, the map T^ £ is ergodic. 

Proof. Assume now that T% £ is not ergodic. Then there exists a bounded 
/ 7^ const, such that T% E f = f, and therefore 

f a M + T a J)=f a J + f 

Since T Qj£ is ergodic, this implies that T a£ f + / = C, or / — C/2 = 
— {T a ^ £ f — C/2). If g = f — C/2, then g is not identically zero, and T a ,e9 = 
—g. Therefore A = —1 is an eigenvalue of T a ^ e and so T Qi£ is not weakly 
mixing. □ 

Proposition 8. Let T be an ergodic transformation on (X,fi), fi(X) = 1, 
and let f G L 1 (X, /i), / fd/i = and S n (f, x) = f{x)+f{Tx) + . . .+f(T n ~ 1 x) 
be its Birkhoff sums. Then either S n (f,x) is unbounded from below for 
almost every x G X, or f is a co-boundary, i.e. there exists a measurable 
g{x), such that f(x) = g(x) — g(Tx). 

Proof. Since T is ergodic, the set of points x for which S n (f,x) is bounded 
from below has measure either equal to zero or one. In the first case, Propo- 
sition is proved, so assume that it has measure one. Then the function 
g(x) = inf S n (f,x) is finite almost everywhere. 

n>l 

We have g(Tx) + f(x) = inf S n (f, x), so h(x) := g{Tx) — g{x) + f(x) > 0. 

n>2 

It is enough to show that f hdji = 0. If g(x) G L 1 (A, /i), then f hdjj, = 
by the definition of h(x) above. If not, then by Birkhoff ergodic theorem, 
for /x-a.e. x G X, 
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lim ^—^ = I hdjj, 



X 



where the integral can be equal to infinity. 
We write 



S n (h,x) _ S n {f,x) g{T n x) - g(x) 



n n n 

Since g(x) is finite almost everywhere, we can choose a set Y G X, such 
that fi{Y) > an d f° r ever Y V ^ Y) \9{y)\ < M for some constant M. Then 
by ergodicity of T, there exists a subsequence n^, such that T nk x G Y, and 
therefore, by Birkhoff ergodic theorem for \x- almost all x G X we have 

lim S nk (h,x) = Hm S nk (f,x) + g{F*x) - gjx) = Q 

which implies J hd[i = 0. □ 

x 

Proof of Theorem^ For any e > choose an a G A(e), so that the map 
T 2 £ is ergodic. Let x G J E and /(x) = n*(x) — (Tq^x) . 
Then the Birkhoff sums for T 2 £ and /(x) are 

x) = f{x) + f(fl £ x) + f(T&) + ... + f(Tl m £ x) = 

= n]{x) - n 2 £ {x) + ...-n 2m (x) 

By Proposition [8j for Lebesgue almost every x G I £ , either there exists 
m G N, such that S mo (f,x) < (and therefore Q e (x,a) < oo) or /(x) is 
a co-boundary. But in the second case, S m (f,x) = g(x) — g(T 2r ^ +2 x) for 
a measurable g(x). Either g(x) < esssup g(x), or g(x) = esssup g(x) on a 
positive Lebesgue measure set. In either case, the ergodicity of T 2 £ , implies 
that for Lebesgue a.e. x, there exists m G N, such that S mo (f,x) < and 
so Q e (x, a) < oo. 

Now let x G 5 1 \J e . Since a ^ Q, there exists n > 0, such that R~ n °x G J e . 
Then for Lebesgue a.e. x E S 1 



QAx.a) 



□ 



3. Limiting distributions 
We now prove theorems [3] and |4| 
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3.1. Notations and the formulation of the main limiting distribu- 
tion result. Let Fe n \ti, . . . t n ) — ¥{em\ > ti, em 2 e > t 2 , ... em" > t n } be 
the joint distribution function of the vector em™ = (em], em 2 , . . . em™) T . 
It is also convenient to introduce 

K{x, a, T) = #{A; £ Z n (0, e^T] : ka + x C I e + Z}, 

the number of times the particle hits vertical walls during the time e~ x T. 
Note that 
(3.3) 

F{em h £ (x,a) > t k , k = 1, . . . , n} = F{J\f e (x,a,t k ) < k - 1 , k = 1, . . . , n}. 

Let xi denote the characteristic function of the interval / C R and 
ipT{x,y) = X(o,i]( x )X[-T/2,T/2](y) be the characteristic function of a corre- 
sponding rectangle. 

Then 

[e-iT] 

Af £ (x, a,T) = ^ Y Xl - ( am + n + %) = 

m=l n&L 



~- Yl x (°^ (^Tt) X[-T/2,T/2]((e x T{am + n + x)) = 

(m,n)el, 2 

= £ «r ((".» + *) (J 



(m,n)eZ 2 

Therefore, 



(3.4) 

A4( a ;,a,T) = #|(m,n)£Z 2 :(m,n + x)^ J i) ^) G ^( T )}' 

where ft(T) = (0,T] x [-1/2, 1/2]. 

Let ASL(2,R) = SL(2,R) x R 2 be the semidirect product group with 
multiplication law 

(M,v)(M',v') = (MI'.vM' + v'). 
The action of an element (M, v) of this group on IR 2 is defined by 



(3.5) w i— > wM + v 

Each affine lattice of covolume one in R 2 can then be represented as Z 2 g 
for some g £ ASX(2,R), and the space of affine lattices is represented by 
X = ASL(2, Z)\ASL(2, R), where ASX(2,Z) = SL{2,Z) K Z 2 . Denote by 
v the Haar probability measure on X. 
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Theorem 9. As e — > 0, the limit of (3.3) exists and is equal to 



(3.6) F<- n \h, ...t n ) = u({g G X : #{Z 2 g nll(t k )} < k-1 (k = 1, . . . , n)}), 
which is a C 1 function IR> — > [0, 1]. 

We define the associated limiting probability density 0( n )(ti, . . . ,t n ) by 



oo oo 



F^(t u ...,t n ) = J ... J <f>^(t 1 ,...,t n )dt 1 ,...,dt n 

t\ t n 

3.2. The reduction of Theorem |3] to Theorem [9l Because of the re- 



lation (|2.2[), Theorem |9| implies the convergence in distribution for the se- 



quences {en^} and {e^} (part (1) of Theorem |3J) . 

Indeed, let k G N and Ii, . . . ,I k be a collection of intervals on the real 
line. Let I = h x . . . x I h c R k . Then 

limP{£^ e / ls . . . ,enj G J*} = limP-Tern?; G B/} = 

" " e— >0 

(j)^(t 1 ,...,t k )dt 1 ...dt k 



BY 

and 



limP{e£ G A, . . . ,eC fc G 4} = limP{em{? G BA^ 1 /} = 
<(>M{t 1 ,...,t k )dt 1 ...dt k 

BA- 1 / 

The convergence for the random variable Q E {x, a) also follows from The- 
orem ini 

Indeed, for any k > 1 let \A k be the characteristic function of the set 

A k = {( yi ,...y k ) G M. k : Vl > 0,...,y k ^i > 0,y k < 0}. 
Then for every e>0we have 

Q £ {x,a)=mm{j eZ + : & < 0} - 1 

Therefore, 

F{Qe(x, a) = k}= F{e£ > 0, . . . e^ 1 > 0, ef* < 0} = 



[em k e G BA-^fc} = J 



e ' 
BA-iAfe 

and by Theorem [9] and the Helly-Bray Theorem ([5], p. 183), there exists 
the limit 
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(3.7) G(k) = limP{Q e (x, a) = k} = [ ^(h, . . . ,t k )dt! . . . 

e^O J 



dti 



Notice that the representation (3.7) implies that G(k) < 1 

k=l 



Proposition 10. 

(3.8) 



k=l 



Proof. Let {r] k } be the limiting process for the sequence {en k }. From the 



explicit description (3.6) of the limiting distribution in Theorem |9| we have 
the following description of the process {rj k }. 
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Figure 2. The horizontal ray through (0,0) generates the 
sequence {— y k } as an orbit of an interval exchange map 



Let g G X be an affine lattice which has no points either with the same 
horizontal coordinates, or on the boundary of the semi-infinite tube TZoo = 
[0, +oo) x [—1/2, 1/2]. The set of such lattices has full Haar measure in X. 
Let us enumerate points of g which lie in TZ^ according to their horizontal 
coordinates: if the coordinates of the fc'th lattice point of g in 7^00 are 
{x k , Vk) = (xk(g), Vk{g)) {k = 1, 2, . . .), then x k < x k+1 for any k — 1,2,... 
Notice, that i/-almost every lattice g has infinitely many points in TZoo. 

The sequence of random variables en k = en k (x,a) w.r.t. the probability 
measure P on T 2 converges in distribution to the sequence T] 1 = ^(g) = 
x 1 (g), and rj k = r] k (g) = x k (g) — x k -i(g) for k > 2 w.r.t Haar measure v on 
X. 
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Therefore, in order to prove (3.8), it is enough to show that for v- almost 



every affine lattice g G X, there exists an even k > 0, such that 
(3.9) T] 1 -r] 2 + r] 3 -. . .-r] h = x 1 -(x 2 -x 1 ) + (x 3 -x 2 )-. . .-(x k -x k -i) < 
We will now show that the sequence y k (g) is an orbit of a certain map 



of an interval into itself, reduce (3.9) to a Birkhoff sum over this map and 
treat it in the way as in Section 2. 

First, we describe the map. Consider set X C IR 2 of vertical segments of 
unit length centered at every lattice point of g. We identify each segment in 
X with the base I — [—1/2, 1/2] of the tube 71^ by parallel translation. Let 
7i : X — > I be the projection, which sends a point on some interval through 
a lattice point to the corresponding point in I. 

Consider a unit speed flow in the positive horizontal direction on M 2 . Its 
first return map to X is a well-defined map T = T(g) of X into itself. We 
define the corresponding invertible map T so that tt oT = T o jr. It 

is easy to see, that the map T is an exchange of three intervals. For //-almost 
every lattice g it has combinatorial type (3 2 1). 

For every y G I, we let ip(y) to be the Euclidean distance between y G 
7r _1 (?/) and its image under T. Clearly, this does not depend on the choice 
of y G 7r~%). 

Notice that the sequence {y k } of the vertical coordinates of the lattice 
points of g in Tloo is related to the map T described above: for k G N, y k — 
— T k ~ x {— yi) (see Figure 2). Also for k G_N, we have ip(—yk) = x k+1 — x k . 



Let — yo = T 1 (—yi). Then the sum in (3.9) has the form (recall, k is even) 



x\ - if>(-yi) + ^(-2/2) - ... - ip(-y k -i) < 
< 4>(-y ) - V(-T(-3/ )) + ^(T 2 (-y )) - ... - ^(T k ~\-y )) 



Therefore similarly to Section 2, the alternating sum (3.9) is reduced to a 
Birkhoff sum for the function f(y) = ip(—y) — T(— yj) and the map T 2 . 

Let the lengths of the interval exchange map T be equal to (Ai, A2, 1 — 
Ai — A 2 ). Denote the simplex of possible Aj's by 

A = {(Ai, A 2 ) I Ai > 0, A 2 > 0, Ax + A 2 < 1} c M 2 , 

and the corresponding interval exchange map of combinatorial type (3, 2, 1) 
by T\ 1 \ 2 ■ The following theorem was first proved by Katok and Stepin in 

ID- 

Theorem 11. For Lebesgue almost every pair (Ai, A 2 ) G A, the map Tx lt \ 2 
of the interval I onto itself is weakly-mixing. 



Similarly to the proof of Theorem [TJ Theorem 11 and Proposition [8] imply 



that there exists a full Lebesgue measure subset Ai C A, such that for every 
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(Ai, A2) G Ai, there exists a full Lebesgue measure subset I' = I'(Xi, A2) C I, 
such that for every y G I' there exists k > 0, such that 

iK-jfo) - ^(-T Al , A2 (-y )) + HTi lM (-yo)) - ... - ^(iti(-w)) < 0. 

Let X C X be the set of lattices, for which the construction above gives 
an interval exchange transformation of combinatorial type (3 2 1). Then X is 
open and u(X) = 1. Notice that for any g G X, the map X : g 1— > (Ai, A 2 , ?/o) 
is differentiable and its differential is surjective. Therefore, the preimage 
of any Lebesgue measure zero set under X has Haar measure zero in X. 
Therefore, the set of lattices g G X, such that X(g) G {(Ai, A2, yo) | (Ai, A2) G 
Ai, y G J'(Ai,A2)} has full Haar measure in X and so (3.8) is proved. □ 



Remark 12. The condition (3.8) is equivalent to the tightness of the family 
of distributions {Q e } as e — > 0. Namely, for any S > there exists N = N(5) 
and E\ = £1(8), such that for e < £% 

N 

(3.11) 1 - 5 < ^>{Q £ (x,a) = k} < 1 

k=l 

3.3. Continuous case. 

Proposition 13. For any s > and 5 > 0, t/iere exists Eq > and fceN, 

(3.12) P{(em e fc < s)} < 8 
for all E < Eq. 

Proof. We have 

P{Em k £ <s} = P{7V £ (x, a, s) > k} 

The limit, as e — > 0, exists and, in view of [7j (p. 1131, first equation), is 
bounded by 

< C s k~ 3 

for some constant C,. □ 



We now prove part (1) of Theorem |4j 

For any N G N and intervals 7i, . . . , In C M, 

PK £ (s 1 )6ii,..,£( £ MeM = 

oo 

= ^ P{ee £ (%) e err# < Sj - < em^ +1 (j = 1, . . . , N)} 



Notice that 



s if < s < em\, 

£ ie + ( — l) fc ( s — £^) if £m k £ < s < £m e fc+1 , 
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and 

k 

e=E(-i) <_i (* 

i=l 

Therefore, by Theorem |9j for every fixed k G Z> we have 

limP{ef e (Sj) G Jj, emj* < < em^ 1 (j = 1, . . . , N)} = 

^ k+1 \ti,...,t k+l )dt 1 ...dt k+1 , 
with k = max(k), and the range of integration restricted to the set 



£^0 



(3.13) 



t kj < Sj < tfc 3 +l, 



Bk — {(ti, . . . , t k+ i) : 
XX-l)*" 1 ^ - i + l)tt + (-!)"' (sj - t k] ) G Aj) 



i=i 



Futhermore, 



e J ^ ^ < *i < (J = 1, • • • , iV)} < 

k G 

maa:(k) > i? 



^ PR £ (s,) G Jj, em^ < Sj < em k e ^ (j = 1, . . . , N)} < 



>o 



< P{emf < si}. 



Part (1) of Theorem [4] now follows from Proposition 13 
For the part (2) of Theorem El we have, 



(3.14) 



keN 



Notice that if Q £ (x,a) = k, then the time which the particle spends in 
the tube is equal to 



T £ = TJx, a) = 2Vl + a 2 (nl + nl + . . . + n*) 



and so, 



(3.15) 



¥{eT £ < s, Q £ (x,a) = k} = 
= F{2eVl + a 2 {n\ + n\ + . . . + n k e ) < s, Q £ (x, a) = k} 
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By Theorem 8, for any s > there exists joint limiting distribution of 

P{a < s, em k (x, a) > t k (k — 1, . . . , n)}, 

as e — > 0, and therefore, of ( 3.15[ ) as well. 
On the other hand, 

F{eT £ < s, Q E {x, a)>k}< F{em k £ < s}, 



and so, Proposition 13 and the convergence of (3.15) imply the existence 
of the limit 



H(s) = limP{eT £ < s} 

£^0 



Also, since 



N 



F{eT e <s}> J2 F i £T e < s > Qe(x,a) = k}, 



k=l 



the tightness (3.11) implies that H(s) — > 1 as s — >• oo. 
This finishes the proof of part (2) of Theorem El 



3.4. The proof of Theorem [9j By ( |3.3 ) it is enough to show that for any 
n E N and any n-tuples (t\, . . . ,t n ) G 
exists the limit 



^ , k = {ki, ...,k n ) G Z^ there 



(3.16) 



limP{A4(a;, a, ^) = kj, (j = 1, . . . , n)} 



£^0 



G X : #{Z 2 £ n ^-)} = fei (3 = 1, • • • , n)}) 

and that G^ n ^(ti, . . . , t n ) is a C^function of (ti, . . . , t n ). 

For n = 1 the convergence in (3.16) was first proved by Mazel and Sinai 
([§])■ It was later reproved and generalized by the third author ([B], [7j) 
using different methods. The proof of the convergence in (3.16) follows the 
one in [BJ. The proof of the regularity of the limiting function is similar to 
the one in [8]. 

We reduce the convergence in (3.16) to an equidistribution result for the 
geodesic flow on X. 

Recall, that the action of the geodesic flow on X is given by right action 
of a one-parameter subgroup of X : 

e~*/ 2 
e*/ 2 



(0,0) 



The unstable horocycle of the flow $* on X is then parametrized by the 
subgroup H = {n_(a;,a)} (;Eia ) eT 2 : 



n-(x, a) 



1 a 
1 



(0,x 



For g G X let Fr(g) be equal to the number of lattice points of T?g in the 
rectangle 1Z(T). 
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Then by (O) 



M £ {x,oc,T) = F T (n_(:r,a)$*) 
with t = -21n(e). 

Theorem 14. [6J For any bounded f : ASL(2, Z)\ASL(2, E) -> E, suc/t 
i/iai the discontinuities of f are contained in a set of v-measure zero and any 
Borel probability measure P, absolutely continuous with respect to Lebesgue 
measure on [0, 1) x [0, 1) 
i i 

hm J J /(n_(x,a)$*)dP(x,a) = J fdv 

ASL{2,Z)\ASL(2,R) 



t— >oo 



Let 

1 iiF t .(g) = kj, (j = l,...n), 
otherwise 



D(g) 



Then D(g) satisfies the conditions of Theorem 14 The convergence in 



(3.16) now follows from theorem 14 applied to the function D(g). 

We now prove C 1 regularity of the limiting function G^ n \ti, . . . , t n ). It is 
enough to consider the case when all tj are different. We also assume that 
all kj > 0. The case when some kj = is similar. 

Let X\ = SL(2, Z)\SX(2, E) be the homogeneous space of lattices of 
covolume one and let v\ be the probability Haar measure on X\. For a given 
yel 2 let 

*(y) = {gex : y e z 2 g }, 

where the action of X on IR 2 is given by the formula (3.5). 



There is a natural identification of the sets X(y) and X\ through 

X(y) = {(M,y) : M G X 1 } 

Under this identification the probability Haar measure V\ on X\ induces a 
probability Borel measure u y on X(y). 
We will need the following two results. 

Proposition 15. (Siegel's formula, [TT]j Let f G L^E 2 ), then 
(3.17) J /(kM)^(M) = J f(x)dx 



kez 2 \o 



Proposition 16. flSjJ Let £ C X be any Borel set; then y — > v y {E fl X(y)) 
is a measurable function from E 2 to E. If U C E 2 is any Borel set such that 
£ C U yeU X(y), then 

(3.18) v{£)< j u y (£nX(y))dy 
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Futhermore, i/Vyi 7^ y 2 G Z7 : -X"(yi) fl X(y 2 ) Pi ^ = 0, £ften equality 
holds in (3.18\) 



Notice that Propositions 15 and 16 imply that if there are two different 
indices 1 < i, j < n, such that 

\j(hi, hj) = {geX : \Z 2 g n {K{U) A Ufa + hi))\ > 0}n 



then 

Therefore, 
(3.19) 



H{g EX : \Z 2 g n (TZ(tj) A 7e(tj + hj))] > 0} ^ 0, 
v{Aij(hi, hj)} = o(||h||) as ||ft|| — > 



G^(ti + hi, . . . , , t n + ft n ) - GW(t 1? . . . , t n ) 



= E G (n) (ii, £ 2 • • • , tj-utj + hj, t j+ i, ...,t n )- G^{h, . . . , t n ) + o(||h||) 
i=i 

= EMfl G x : l z ^ n = fc i - h \ z2 9 n ^fo + h i)\ = k 

|z 2 #n7e(t;)| = ki,i^j}- 

-u{geX : |Z 2 ^ n 7^(^-)| = \Z 2 g mZ{tj + hj)\ = kj + 1} 

|z 2 0nft(*i)| = fc,i^j) + o(||h||) 

Consider a single term in the expression above. 
Let 



5j = ^(/ij) = {s6l: \Z 2 g fl 72.(tj-)| = 



|z 2 # n + hj)\ = kj + 1, \Z 2 g n 7^(^)1 = h,i^ j}, 



and let {7 = 7£(ij + ft,,-) \ H(tj). Then by the proposition 16 



1/2 



= y i/ y (£j n X(y))dy = J J V( x , y ){£j n X(x,y))dxdy 



u 



tj -1/2 



Therefore, by proposition 15 



. lim ft T-K^i('li)) 
hj— >0 ftj 



1/2 

y ^eli : |Z 2 <?n(ftfo)-fo,y))| = fci, (i = l,...,n)})dy 



-1/2 



For every fixed y G [—1/2,1/2] continuity of the expression under the 



integral sign with respect to (£i, . . . ,t n ) again follows from Proposition 15 



It is clearly uniform in y and therefore the integral is continuous with respect 
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to (ti, . . . , t n ). Each term in (3.19) can be treated in a similiar way and this 
proves C 1 regularity of the function G^ n \ti, . . . , t n ) and finishes the proof of 
Theorem |9l 

4. Proof of Theorem [2] 



We now deduce (1.1) from part (2) of Theorem [3] 

Consider the unfolding of the tube to IR 2 obtained by the reflections from 
the horizontal boundary of the tube. Let pk = (£* , ( £ ) be the position of 
the particle at the moment of fc'th reflection from the wall in this unfolding. 




Figure 3. An unfolded trajectory. In this example, Q £ 
3, LCJ = 2 and n\ = 2 



, ti* = 1, n 6 £ = 3, rf £ > 4 



Then 
and 



g = n2-n2 + ... + (-i; 



C* = Vin + a{nl + n\ + . . . + n k £ ) 



At the moment of the exit from the tube, the vertical coordinate of the 
particle is 



(4.20) 
Let 



and let 1 1 • 
Then 

Therefore, 



C = 2(y in + an\ + an z e + . . . + anf") - y h 



z — Vin + OLn\ + anl + . . . + an, 
denote the distance to the nearest integer. 

| \yi n + an l £ \ \ < e/2, \ \an l £ \ \ < e for i > 1 



(4.21) 
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Notice that t> ou t = — v- m , if both the number of reflections from vertical 
walls and from horizontal walls is odd. The former is obviously odd at the 
moment of exit. The number of reflections from the horizontal walls is equal 
to the integer part |_CJ • 



If z — [_z\ < 1/2, then by ( |4.20[ ), |_CJ is odd provided that 2||z|| < y- m , and 
if z — [z\ > 1/2, then |_CJ is °dd provided that 1 — 2||z|| > y in . 
By (4.21) this is the case, when 



eQ £ < min{y in ,l - y m } 

By the assumption, the probability measure P on the initial conditions 
(yi n , a) is absolutely continuous with respect to the Lebesgue measure, there- 
fore for any fcGN, 

P{Q £ = k, ek < min{?/ in , 1 - y m }} = 
HQs = k}-P{Q e = k, min{y in , 1 - y in } < ek} -> G(k) as e -> 



Together with the tightness condition (3.11) this implies 
P{sQ £ < min{?/ in , 1 - y in }} -> 1 as e -> 

and so, 

P{^out = -^in} -> 1 as 6 -> 0. 

Note that the existence of the limiting probability distribution for {Q £ } as 
e — > also implies that for any 5 > 0, 

P{|j/out - 2/inl > 5} -> as £ -> 0. 

Indeed, after each reflection from a vertical wall, the particle backtracks 
itself with an error at most e, so at the moment of exit it backtracks the 
incoming trajectory with total error of at most eQ £ . 
This finishes the proof of Theorem [2] 
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